Nonlinear double Compton scattering in the full quantum regime 
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A detailed analysis of the process of two photon emission by an electron scattered from a high-intensity laser 
pulse is presented. The calculations are performed in the framework of strong-field QED and include exactly 
the presence of the laser field, described as a plane wave. We investigate the full quantum regime of interaction, 
where photon recoil plays an essential role in the emission process, and substantially alters the emitted photon 
spectra as compared to those in previously-studied regimes. We provide a semiclassical explanation for such 
differences, based on the possibility of assigning a trajectory to the electron in the laser field before and after 
each quantum photon emission. Our numerical results indicate the feasibility of investigating experimentally 
the full quantum regime of nonlinear double Compton scattering with already available plasma-based electron 
accelerator and laser technology. 
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Electromagnetic radiation by accelerated charged particles 
is one of the most fundamental processes in physics and it 
is exploited experimentally for different purposes, spanning 
from the generation of coherent x-ray fields [1] and of even 
multi-GeV photon beams [2] to medical applications. From a 
more fundamental point of view and by limiting to the case of 
electromagnetic driving fields, the process of electromagnetic 
radiation by accelerated charges has played a crucial role for 
testing the validity of classical electrodynamics and of QED. 
In QED the emission of radiation by a charge, an electron for 
definiteness (mass m and charge e < 0, respectively), is de- 
scribed as the emission of quanta of the electromagnetic field, 
i.e., of photons. As in classical electrodynamics an electron 
emits only if it is accelerated, quantum mechanically pho- 
ton emission can occur only if the electron absorbs at least 
one photon. (Linear) Compton scattering, i.e., the emission 
of a single photon by an electron via the absorption of an- 
other photon, has been investigated theoretically and experi- 
mentally since the formulation of QED itself, and it has pro- 
vided relevant information on the structure of the theory at 
high energies [3]. The related process in which two photons 
are emitted ((linear) double Compton scattering) has also been 
investigated both theoretically [4] and experimentally [5]. 

When an electron is driven by an intense electromagnetic 
field, the emission of photons may occur with the absorp- 
tion of many photons from the field. High-power lasers are 
sources of strong electromagnetic fields with unprecedented 
intensities [6] and they represent an irreplaceable tool to test 
the high-intensity or strong-field sector of QED, as comple- 
mentary to the high-energy one. In strong-field QED the back- 
ground field is so intense that it has to be taken into account 
exactly in the calculations. For a laser field approximated as 
a plane wave, this is required when the Lorentz- and gauge- 
invariant parameter ^ = \e\£L/iTicujL is of order of or larger 
than unity, where £l is the electric field amplitude of the wave 
and ojl its central angular frequency. The threshold £, ^ I cor- 
responds to an optical (hojL ~ 1 eV) laser intensity of about 
10^* W/cm^, which is four orders of magnitude smaller than 
available optical laser intensities [7]. The emission of a single 



photon by an electron in the field of a plane wave (nonlin- 
ear single Compton scattering (NSCS)) has been investigated 
theoretically since the invention of the laser (see the recent re- 
view [6]) and a pioneering experiment performed at SLAC [8] 
has confirmed the predictions of QED at ^ < 1. More recent 
studies on NSCS have been focused especially on the high- 
intensity regime ^ ^ 1, where a large number (~ of laser 
photons is absorbed by the electron during the emission pro- 
cess, and on finite-pulse effects [9-12]. NSCS is found to be 
characterized also by the so-called quantum nonlinearity pa- 
rameter x = {{kLp)/mojL){£L/£a)- Here, = {ujl/c, Ul) 
is the four-wavevector of the laser photons (|fci| = lol/c), 
— {e/c,p) is the initial four-momentum of the electron 
and £cr — m^c^ /h\e\ = 1.3 x 10^^ V/cm is the critical field 
of QED [6] . The parameter x controls quantum effects like the 
photon recoil and at x ^ 1 the NSCS spectra coincide with 
the classical ones [13]. The emission of two photons by an 
electron in the field of a plane wave (nonlinear double Comp- 
ton scattering (NDCS)) has also been studied in the Uterature, 
with an emphasis on the correlation of the two emitted pho- 
tons [14] and on the relative yield between NSCS and NDCS 
[15]. Both these studies have been focused on the radiation 
regime where ^ ^ 1 such that the electron absorbs only a few 
photons from the laser field ("quasilinear" regime), and where 
X <C 1, such that quantum photon recoil was negligible. 

In the present Letter we investigate NDCS in the full non- 
perturbative quantum regime where ^ 1 and x ^ 1- By 
including exactly the effects of the plane wave, we show that 
the emission spectra in the full quantum regime substantially 
differ from those at x ^ 1 as a result of the photon recoil. 
We explain the new features of the emission spectra, by de- 
veloping a quasiclassical approach, based on the possibility of 
assigning a trajectory to the electron in the laser field, with dis- 
continuities in the electron energy due to the quantum emis- 
sion of photons. Finally, by means of a numerical simulation, 
we show that the full quantum regime can be entered already 
with experimentally demonstrated laser intensities of the order 
of 10^^ W/cm^ [7] together with presently available plasma- 
based electron accelerator technology [16]. 
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FIG. 1. Tree-level Feynman diagrams of NDCS in the Furry picture. 
The double solid lines represent Volkov states and propagators. 



We consider a linearly-polarized plane-wave field, de- 
scribed by the four-potential ^^(j]) = {£l/^l)^i,'^a{v)^ 
where is the wave's polarization four- vector, and ipAiv) 
gives the field's temporal shape, which depends on the space- 
time coordinates a;^ only via the invariant phase ij = 
(units with h = c = 1 are used throughout). For a generic 
four- vector a^^ = {a", a) it is convenient to introduce the 
light-cone representation a'' = (a"*", a~, a-*-), where = 
(a°±all) /V^, with a" = ■ a/ujL, and where = 
a — o'^kL/ujL- In the Furry picture of QED, the electron 
wave functions are solutions of the Dirac equation in the pres- 
ence of the considered background field, i.e., Volkov wave 
functions in the case of a plane-wave field [17]. The Volkov 
wave function for an electron with four-momentum and 
spin quantum number a outside the plane wave has the form 



'^p,a{x) = Ep{x)Up^„/y/2e, where Up^cr is a free bi-spinor, 
where a unity quantization volume is assumed, and where the 
Ritus matrices 



Ep{x) 



1 + 
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7''ap for the 



have been introduced, with the notation ( 
Dirac matrices 7'^. We consider an electron with initial (fi- 
nal) four-momentum = {ei,pi) (pj = {ej^pj)) and spin 
quantum number ai (cr/), which emits two photons with four- 
momenta fcf and and with polarization four- vectors e^^ 



and respectively. The scattering matrix element Sfi of 

^fi ' "^/i ' where (see Fig. 1) 



-k2,\2 

this process can be written as s9'^ '^^^^ 



(2) 



and is obtained from S'j^^ by exchanging the photon in- 
dices 1 and 2. In Eq. (2) we have introduced 'i'p^„{x) = 
*p,cr(^)7°' the dressed electron propagator [18] 



G{y,x) 



lim 

e->0 



^ Ep{y) J . ^. Ep{x) (3) 
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with Ep{x) = 7°£^(a;)7°. It is clear that it is sufficient to 

analyze here only the quantity s'j^K The structure of the Ri- 
tus matrices allows to perform the integrals with respect to 
the coordinates a;+ and x-^ (y+ and y-^), which provide the 
energy-momentum conservation laws = k^ + and 



ki + P i p-^ = k2 + pf and p 



Pf) at 



vertex x (y). One of these sets of conservation laws can be 
employed to perform three integrals in the electron propaga- 
tor. The only remaining integral I{y~ , a;~ ) in Eq. (3) is 



I{y-,x-) = lim / ^ ^+7^ e-^^(---^ 
' e^oj 2n P+-P+ +ie 

= fiLS{x~ - y~) -{{■^t+m)Q{y~ - x~), (4) 

where = k'^/uj^, is the step function, and where we 
have introduced the "transitional" four-momentum with 
light-cone coordinates p^ — p~ — k^ ~ + pj , p^ = 

pj- — ki = k^ +pf andpt^ = (p^^ + m^)/2p^. Note that 
the four-momentum p^ is on-shell, i.e., it fulfills the condition 
p1 = m? . The above decomposition allows to write the quan- 



tity Sf^ as S 



:(i) 

fi 



{2n)H{pT - k^ 



k^ - pJ)S(^\pI - 

ki -k^ - pf) Y.l,s=a((^rfrSr,s + K,sfr,s), whcrc the co- 
efficients ar and fe^.s are matrix factors given by products of 
a large number of Dirac matrices, whose exact form is not 
needed here. In fact, all the dynamical information on the 
process is contained in the functions 



U=j dj^VAiv) exp{-i [SM + Sy(v)]}, (5a) 
fr,s = j dr]^dr]y@{r]y - r?^) i'AiVxWAiVy) 



X exp{-i [S'x(77x) + Sy {%)]}, 



(5b) 



where Sj^/y (77) = Jq dr]'[ tpA{v')+l3x/y'^Ai'n')+7x/y], 
with ax = -m^lipreL) / {kLPi) - {pteL) / (kLPt)], Px = 
-m^^'^{kLki)/2{kLPt){kLPi), 7x = -{km)/{kLPt) and 
with ay, j3y and obtained from a^, Px and 7^, respec- 
tively, with the substitutions p'^ pp p^ pf and fcf 
k2. The integrals, which do not contain the shape function 
'^Aiv) in prefactor, diverge. These divergences can be 
avoided by employing the identities [10, 15, 19] 7y/oj = 
-ifj - ay fij + (5y f2,j, 7x fi,o = i/i - fi,i + Px fi,2 and 
ilx + 7y)/o = -{ax + ay)h- {Px + Py)f2. The differen- 
tial average energy emitted dE, summed (averaged) over all 
outgoing (incoming) discrete quantum numbers is given by 



dE 



LOi +UJ2 Cp'Pf TT Cpk. 



n 



?(2) 
'/i 



(6) 



where {<t, A} = ct^ , ct/ , Ai , A2 . Note that the three 5-f unctions 
contained in (and equally in S^^^) can be exploited to 
perform the integrals in p/. 

We will consider the electron-laser collision in a reference 
frame where the laser pulse propagates along the positive z- 
axis (n^ = (1, 0, 0, 1)), it is polarized along the cc-direction 
(e^ = (0, 1, 0, 0)) and the electron is initially counterpropa- 
gating with respect to the laser beam, i.e., pi = (0,0, --(iiEi), 
with the electron's initial velocity — /?« < 0. In order to 
stress differences to existing treatments valid for monochro- 
matic fields [14], the laser's temporal pulse shape is modeled 
by the function iPa{i) = sin'* (77/4) sin(r7) for 77 e [0,47r] 
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and zero elsewhere, corresponding to a two-cycle pulse of ap- 
proximately 5-fs duration at ujl ~ 1-55 eV. We first consider 
a laser system with peak intensity = 5 x 10^^ W/cm^, 
corresponding to a classical nonlinearity parameter w 15. 
In this regime the process exhibits a highly nonlinear depen- 
dence on the laser amplitude, as about f ^ 3000 photons 
are estimated to be absorbed by the electron from the laser 
field [10, 18]. Moreover, in order to highlight the qualitative 
differences between the full quantum regime x ~ 1 consid- 
ered here and the already-studied regime at x ^ 1, we first 
consider a numerical example within the latter Thus, we set 
£i = 40 MeV, which corresponds to x = 5 x lO"'^, ensur- 
ing that photon-recoil effects are negligible. We choose to 
observe one photon at 6'i = tt — 9q/2, with 6*0 = m£^/ei, 
and the other one at the two different polar angles 62 = Oi 
(see Fig. 2a) and 62 = — 1-1 ^0 (see Fig. 2b). Also, we 
choose (pi = Tr,ip2 — in both cases as azimuthal obser- 
vation angles. We observe that for our pulse shape ipAiv)' 
1) the emission cone of NSCS is determined by the condi- 
tion TT — 9 < 0.8 6*0 on the polar angle 6 (the pref actor of 6*0 
is given by the numerical maximum of the function ipAiv)) 
[20], and 2) the NSCS emission spectra at {9, = 0) and 
at {9, (j) = tt) coincide. Therefore, if 92 = 9i both photons 
are observed within the emission cone of NSCS, whereas if 
02 = TT — 1.1 9q one of the photons is observed outside this 
cone. Figures 2a and 2b show that there is emission of radi- 
ation only within the NSCS emission cone. This feature can 
be quantitatively understood by virtue of a stationary-phase 
analysis [10, 18]. In fact, in the regime m^,£i ^ m the 
phases in the integrands in Eqs. (5) are of the order of and 
thus very large [10, 18], and the saddle-point method can be 
applied to evaluate the functions fr and Jr^s- In the case of 
NSCS, the stationary points for a given observation direction 
correspond to those phase instants where the classical elec- 
tron velocity points towards that observation direction [10]. 
Accordingly, the overall emission cone corresponds to the an- 
gular region that is spanned by the classical electron's velocity 
vector along the complete classical electron's trajectory in the 
laser pulse. Thus, the electron propagation in the laser field is 
quasiclassical and the photon recoil is the main quantum effect 
to be accounted for [21]. Now, we have checked numerically 
that the emission spectrum is dominated by the contribution 
proportional to the functions fr^s- These integrals, accord- 
ing to the stationary phase method, can be approximated as 
fr.s ~ Ei,nQ(%,» - Vxd)fi!{Vv,n)fsi'nx,i)- Here, the in- 
dices / and n run over all stationary points, which are found 
as solutions of the equations dS ^ / y{r]) / drj\,-i^f^^ ^ =0 and 

/."/^ = / dT^VAi-n) exp[-i5,/,(77)] (see Eq. (5b)). The 
above-approximated expression of fr^s can be interpreted as 
a two-step emission process in which the electron first emits 
a photon with four-momentum fc^ changing its own four- 
momentum from the initial one to the transitional one p^, 
and then it emits a photon with four-momentum fcj chang- 
ing its own four-momentum from the transitional one to 
the final one (recall that the total amplitude also contains 
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FIG. 2. (color online) Two-photon energy emission spectra 
d£/nf^idwidf74eV"^ sr"^] atx = 5x10"-^ observed at = tt- 
6*0/2, and at 62 = O1 (part a)) and at 6I2 = tt - 1.1 So (part b)), with 
^0 = 0.19 rad. The other numerical parameters are given in the text. 
Part c): NSCS emission probability dW^^^^ /Aujdn\€W^^ sr"^] at 
{9\,(j)i). Part d): classical electron trajectories (the instantaneous 
electron energy is color-encoded) with initial four-momenta and 
p^, joined at f]x,\. Since as a typical photon emission energy the 
value 10^ eV has been chosen (see parts a) and c)), then ~ 
and the two trajectories are indistinguishable. Emax is the maximum 
electron energy in the laser field. The crosses mark the points of the 
trajectory, where the classical electron's velocity is along {6\,<j)\). 



a term with the photon indices 1 and 2 exchanged). In ad- 
dition, one can picture this dynamics as a succession of two 
classical trajectories with initial four-momenta and p^, re- 
spectively, continuously joined at a point corresponding to the 
phase f]x,i where the photon with four-momentum /cf is emit- 
ted. In Fig. 2d we show in principle a pair of such classical 
electron trajectories joined at the point marked with a cross 
and labeled as fjx,i where the electron propagates along the 
observation direction (6*1, 0i). However, since in the regime 
X ^ 1 the recoil is negligible (see Figs. 2a and 2c and recall 
that £i — 40 MeV), the two four- momenta andpf are prac- 
tically identical and the two trajectories are indistinguishable 
(in Fig. 2d we subtracted the momentum of a typical photon 
of energy wi = 10^ eV emitted towards (0i,0i) from the 
initial electron momentum, see Figs. 2a and 2c). The same 
occurs if the two trajectories are joined at the other saddle 
point f]x,2 (see Fig. 2d), where the electron's velocity is again 
along the observation direction (6*1, 0i). From an analytical 
point of view, by inserting the parameters and (i^ (ay and 
(3y), the stationary-point equation at the vertex x (y) evaluates 
to the approximate conditions ipA (Vx) = ^1/2 (tpAiVv) — 
M2 - (wi/£.)(l/2 + M2), with M2 = - 92)/9o). If 
photon recoil is negligible (cji ^ Si), we observe that the 
condition for fjy has a real solution only if Az92 < 0.8, i.e., 
only if 02 lies within the NSCS emission cone. This explains 




the absence of emission outside this cone in Fig. 2b. Also, 
by comparing now the frequency distribution of the NDCS 
emission spectrum in Fig. 2a with the NSCS emission prob- 
abihty dW^^'^^ /dujdQ [9-11] in Fig. 2c, it is apparent that 
the NDCS spectrum for x ^ 1 corresponds to an emission 
probability given by the "product" of two independent NSCS 
probability distributions for each photon (note that the max- 
imum in the probability in Fig. 2c is at lower energies than 
in the two-photon emission spectrum in Fig. 2a, as the latter 
contains an additional factor uji + 0^2). Since the formation 
length of NSCS at ^ ^ 1 is much smaller than the laser's cen- 
tral wavelength = 2tt/ujl [6], if recoil effects can be ne- 
glected, the two-photon emission process can be described by 
two independent NSCS events [22]. Accordingly, we checked 
that if VKNDcs (i^nscS) jj^g j^j^j ^j^q^ (NSCS) emission 
probabihty, then W^^dcs ^ (•|^nscS)2/2 [22]. 

The physical situation, however, changes substantially if we 
enter the full quantum regime at x « 1. In order to investigate 
this regime, we set Si — 2.5 GeV, as was deemed achiev- 
able by current technology [23], and — 3 x lO^^W/cm^ 
(^ « 37), resulting in x = 1.1, and keep all other parameters 
unchanged with respect to the above example. By observing 
the two emitted photons at the same emission angles as before 
(see Figs. 3a and 3b), we note that: 1) the quantum mechani- 
cal cutoff energy for the sum of the emitted photons' energies, 
approximately given by the equation Ldi+uj2 = Ei, is well ap- 
proached; 2) since in the quantum regime the electron loses a 
substantial part of its energy after the emission of the first pho- 
ton, the asymmetry in the energies of the two emitted photons 
(see Fig. 3a) is much more pronounced than at x ^ 1 (see 
Fig. 2a); 3) the electron also emits outside of the NSCS emis- 
sion cone (see Fig. 3b and note that Oq — 7.6 x 10^'' rad for 
the present numerical parameters). This last feature is particu- 
larly important as it can be exploited to measure NDCS in the 
full quantum regime. In order to explain it qualitatively, we 
show in Fig. 3c the classical trajectories for initial momenta 
Pi (solid line) and pt (dashed line), obtained for the emission 
of a photon of energy uji = 0.8 GeV and momentum along 
(^1, (pi) atfjx^i (the analysis in the case in which the electron 
emits at fjx,2 is analogous). At this value of uji the emission 
probability is maximal. The derivative of the two trajectories 
is continuous at fj^^i (see Fig. 3c), as the photon is assumed to 
be emitted with momentum parallel to the instantaneous elec- 
tron's velocity. However, the color-coding shows that the elec- 
tron energy discontinuously decreases at the emission point 
due to photon recoil. The abrupt decrease in energy induces a 
stronger deflection of the electron trajectory in the laser field 
after the photon emission. Thus, the resulting emission cone's 
opening angle increases and the emission outside of the NSCS 
emission cone becomes possible. The same conclusion can 
be drawn analytically, based on evaluating the saddle-point 
equation ipAiVy) = Ai?2 - (wi/ei)(l/2 + Ai92) < 0.8, with 
At?2 = {tt — O2) / 9o- In fact, solving the latter inequality for 
62, we obtain an expression for the cutoff angles of the second 
photon emission as a function of loi (see Fig. 3d, with 02 = 
mapped to 6*2 > vr and 4)2 — mapped to 62 < tt). The 




FIG. 3. (color online) Two-photon energy emission spectra 
dE /IQ^^A^JJ^AQ.^[e^^-^ sr"^] atx ~ 1.1 observed at Si = 7r-eo/2, 
and at 62 = 9i (part a)) and at ^2 = tt — 1.1 So (part b)), with 
6() — 7.6 X 10"'' rad. The other numerical parameters are given in 
the text. The solid white lines correspond to the cutoff-energy equa- 
tion uj\ -\- UJ2 = Ei. Part c) the two classical electron trajectories with 
initial electron momentum pi (solid line) and pt (dashed line). The 
color-encoded line shows the actual electron trajectory for a photon 
with energy cji = 0.8 GeV and momentum along {9i , 01 ) emitted at 
f]x.i. Part d) Emission opening angle for the second emitted photon 
as a function of uji (light shaded area) compared with the emission 
cone for NSCS with initial electron momentum pi (dark stripes). The 
vertical line indicates the value toi = cut described in the text. 



photon-energy threshold luI, beyond which emission along 
62 = TT + 1.1 6*0 becomes possible (see Fig. 3b) is well re- 
produced by this analytical prediction of the angular emission 
region (see Fig. 3d). In this way, measuring MeV photons 
outside of the NSCS emission cone would reveal a NDCS sig- 
nal, where NSCS is exponentially suppressed and negligible. 
We have estimated numerically that the emission probability 
outside of the NSCS emission cone is of the order of 0. 1, indi- 
cating the observability in principle of the process, employing 
laser-generated electron beams, which typically contain of the 
order of 10^ electrons [16]. 

In summary, we have analyzed for the first time nonlinear 
double Compton scattering in a plane-wave field in the full 
non-perturbative quantum regime, where a large number of 
laser photons are absorbed by the electron during the emission 
process, and where photon recoil substantially alters the elec- 
tron's dynamics. This alteration is, in particular, responsible 
for the emission of photons along directions, where nonlin- 
ear single Compton scattering is suppressed, which provides 
a simple way of detecting double Compton scattering in the 
full quantum regime already with available laser and electron 
acceleration technology. 

The authors acknowledge useful discussions with A. Ilder- 
ton, D. Seipt and B. Kampfer 
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